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Abstract :FoUowing Dahl’s method an exact Runge-Lenz vector M with two components 
Ml and M 2 is obtained as a constant of motion for a two particle system with charges ei 
and e 2 whose electromagnetic interaction is based on Chern-Simons electrodynamics. 
The Poisson bracket {Mi,M2}9^Lz but is modified by the appearance of the product eie 2 as 
central charges. 
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A characteristic feature of the non-relativistic Kepler problem is that there exists 
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apart from the conservation of energy(E) and the angular momentum vector (L) another 
conserved quantity ,namely the Runge-Lenz vector (A). A textbook^ derivation of A 


begins with the equation of motion of a mass m under the central force F = f( r)—, 

r 

dp r 

namely, — = f(r)— , leading to 
dt r 


d 


d 


Ot / \ 2 ^ 

-{pxL]=-,nf(r)r - 


frl 

yr J 


( 1 ) 


dL 2 

as — = 0. With r f (r) = a constant(say - k),implying thereby an inverse square law of 
dt 


r 

force, eq. (1) immediately yieldsA = pxL - mk — ; this is the Runge -Lenz vector. With 

r 

the definition K = A, one easily obtains the Poisson bracket relations 

{K.E] = 0.{L.E) = 0.{L,.L] = t„,L,.{L,.K] = e,„ K,.{k,.k] = e„. L. (2) 

Thus there exists an internal symmetry associated with the non-relativistic Kepler problem 
with the invariance group being isomorphic to the 4-dimensional rotation group O 4 . 
Until recently,the presence of this internal symmetry had not been tied to a generally 
accepted invariance principle. In other words,the phenomenological derivation imphed in 
eq. (1) begs the question of whether the Runge-Lenz vector has a deeper physical origin. 
Specifically, the question is if there is a space-time transformation,the invariance of the 
Lagrangian for the Kepler problem under which directly leads to the conservation of A. 

An affirmative answer to this question was recently obtained by Dahl ^ by regarding the 
Kepler problem as the zero-order description of a relativistic two-body problem;or,as 
emphasized by Dahl it is absolutely necessary to investigate the relativistic two-body 
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problem in order to discover the connection between the dynamical symmetry of the non- 
relativistic Kepler problem and special relativity. 

Of special relevance to this paper is the fact ^ that for the Kepler problem the angular 
momentum vector L = r x p is conserved so that r . L = 0. Thus r always hes in a 
plane whose normal is parallel to L. Since the motion of the particle is planar because 
of symmetry considerations it is appropriate to ask: Suppose the motion of the particle 
was de facto planar and not due to symmetry considerations as in the case of motion under 
a central force,could one stiU obtain a Runge-Lenz vector that is conserved in this 2+1 
dimensional case? Happily,we derive in this paper following Dahl’s method ^ a bonafide 
Runge-Lenz vector (see eq. (23) below)associated with the relativistic Lagrangian with 
Chern-Simons interactions ' for a two-particle system consisting of masses mi and m 2 
with charges ei and e 2 respectively, given by 


L = -, 


mjc 


1 - 


.2 

Xi 


\l /2 


■m2C 


1 - 


.2 

X2 


\l /2 


6162 k.rxr 


Tic 


( 3 ) 


with r = xi - X 2 ,xi (and Xj) and X 2 (and jC 2 )being the position (and velocity) vectors of the 
masses mi and m 2 and k a fictitious unit vector orthogonal to the plane of xi and X 2 . 

It is easy to obtain (3) starting with the electromagnetic potentials (x,t) ,p = 0,1,2 for 


Chern-Simons electrodynamics in the radiation gauge with 'V.A = 0 ,namely. 


(x,t) 


e (jc 
27tc 


-R(t))xR(t).k 

x-R(t)^ 


(4a) 
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e (x-R(t))xk 

A(x,t) = —- -(4b) 

2^ \x-R(t)\ 

the charge e being located at R (t) .Needless to say so, eqs.(4) are the 2 + 1 dimensional 
counterparts of the corresponding 3 + 1 dimensional (x ,t ) given by eqs.(26.19) and 
(26.20) in Fock ^ for example. 


Under an infinitesmal Lorentz transformation given by Siv = -vt + —^x ,it 

c 

is simple to check that the change in the Lagrangian L which is defined by 

dL ^ dL ^ dL dL 

OL = - -.QV; +T^- .0X7 +i:—.OX] 

oX j ^*^2 ^X j ^*^2 


d i ^ 1^2 k.drxr 

— -mjXj.vyj-m2X2.vy2+— - 

dt[ nc \r\ 


without the use of the equations of motion for the two masses mi and m 2 .Here yi and yi 


are given by 1 ■ 


and 1 ■ 


respectively. Since 5L works out to a total 


differential in (5) without the use of the equations of motion, it is clear that the action 
S = \dt L is unaffected by the transformation from jc ^ jc + dv thus making the 
Lagrangian (3) Lorentz invariant. 

The constant of motion is now obtained by using the equations of motion for mi 
and m 2 to rewrite the first equahty in (5) as 


ddL 

oL = — .OK j + .OK 2 

dtyoXj 0X2 


From the second equahty in (5) and (6) one thus obtains 
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dt 


^ dL dL eje2 k.brxr^ 

— .bXi+—.bX2 +miXi.VJi+m2X2.Vj2 -^ 

axj 0X2 T^c r 


= 0 (7) 


dL 

With the canonical momenta defined by pi = it is easy to obtain from (3) 

OXi 

_ 1 _ 1 I 12 

Pi =mjX]yj +g and P 2 =m 2 X 2 y 2 ~8 ,where 7Cc|r| g = eje 2 (r xk ) and thus rework 

d , , 

(7) as — [K.v) = 0 with the constant of motion K given by 

dt 


K = -t P + m ]X ly i^+m 2 X 2 y 2 ^ 


( 8 ) 


the total momentum P being defined as P = pi + P 2 -Eq. (8) is the counterpart of eq. 
(20) in Ref. (2);also while the latter is derived from the Darwin Lagrangian ^ for the 
electromagnetic N-body problem and is correct in the 1/c^ approximation , eq. (8) above 
is exact. Note that K depends exphcitly on time except in the centre-of -momentum 
system which is the Lorentz frame where P = 0. 

For the Lagrangian (3) the Hamiltonian H and the angular momentum vector L are 
given by 

H = nij y~j'+ m 2 y ~2 , L = Xj XPi + X 2 XP 2 (9a) 

and in the P = 0 frame they work out to 

H=mic^ fi+m2C^ f2, L = rxp = k(r,p2 - r2Pi ) = kL^ (9b) 

with Pi and p 2 being the x- and y-components of the relative momentum vector p 

m 2 Pi - mi P2 

defined by p =-. 

m; +m2 


In deriving the Hamiltonian in (9b) we go through the following steps: 
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—1 —2 ^ \ 1 ^ 
i)Since m; jCyY; ^, therefore Y; = 1+ ^ 2 \Pi ~S\ . One can now express pi 

mj c 

and p 2 in terms of p and P via the relations 


Pi = 


nij 


nil +^2 


P + p , 


P2 = 


m2 


mi +m2 


( 10 ) 


and find that in the centre - of - momentum system with P = 0 




B 


2 A 


i + ' 2 2 

V mi c J 


(11a) 


with B = p - g . 


ii)The above calculation can be repeated with m 2 x 2 Y 2 ^ = P 2 + S to obtain in the case 


when P = 0 the result 


12 = 


1 + - 


B 


2 \ 


2 2 
m2 c 


(11b) 


We shall re-label y 1 and y 2 '^ in eqs.(ll) as f 1 and f 2 respectively below.Using these 
definitions one can now rework K in eq. (8) in the centre-of-momentum system with P 
= 0 as 


K = (mifi +m 2 f 2 )R + \i(fi-f 2 )r 


( 12 ) 


with the reduced mass |i and the centre - of - mass position vector R defined by 


mi m 2 miXi+m2X2 

|i =-and R = -. 

mj +m2 mi + m 2 

As emphasized by Dahl ^ the c. m. vector R should now be eliminated from eq. (12) 
to make K a proper dynamical function and as in Ref. 2 we shall do this with the help of 
the the equation of motion for R in the P = 0 hmit,namely, 
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dR (dH 
dt [ dP 


il dH ii dH 

— _!-- 

P=o [rni^Pi rrn ^Pi 


Pl-P’P2 ~~P 


One should now recast the Hamiltonian H in (9a) in terms of the canonical momenta pi 
and p 2 and it is given by 


f I ,2 ^1/2 ( I ,2 

2 2 \^l\ P2I 

H = mjc g]+m 2 C g 2 , gj = 1 + ’^2 = ^ + 


2 2 
nij c 


2 2 
m2C 


with Si = Pi - g and S 2 = P 2 + g • On calculating the partial derivatives required in eq. 


(13) it is easy to get 


(fni+tn2] — = (yi-y2)B 


From the Hamiltonian (9b),the equations of motion are: 


TB = r 


(16a) 


Yi Y 2 

with T = — + — , Y 1 and y 2 being defined by eqs.(l !)• Since B = p - g it is clear that 
ni] m 2 


k xr = T{k X p + hr) 


(16b) 


where Tier h=-eje 2 . With eq. (16b) the other equation of motion ,namely, 


dH I r.kxp 

T— =Th k X p-r h + 2 —^— 

dr r 


(16c) 


yields 


p = h kvr-2--ykvr.r 


Eq. (17) easily leads to 
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d ( \ dB 

—\p + hrxk) = —— = 0 
dt dt 


( 18 ) 


Thus B is a constant in time and so are f i and f 2 in eqs.(ll). Returning to eq. (15) we 


get 


dt 


^ Y 2 -Y 7 ^ 

R + ———B 

mj+m2 j 


0 


(19) 


on using eq. (18). Thus the c. m. position vector R(t) is given by 

(mi+m2)R(t) = B(yi-y2)t + Ro 

where Ro is a constant. Eq. (20) is the counterpart of eq. (43) in Ref. 2;however we 
must emphasize here that eq. (20) above is exact in that it has been derived using the 
equations of motion given by (16a) and (16c) without recourse to approximation. In 
contrast,the solution for R(t) given by eq. (43) in Dahl’s paper ^ has been obtained in 
the 1/c approximation to the equation of motion as given by eq. (39) therein. On using 
(20) in (12) one obtains the exact form of the vector K as 


K 


^ mi-m2+m2 f2 Ji-mj fj Y 2 ^ 
nil + ^2 


mjfj+m2f2 

\tB + ii(fi-f2)r+R, --- ( 21 ) 

ni] +m2 


With the Hamiltonian H in (9b),eq.(21) can be reworked as 


K 


H nil —ni’y 

Rn+— - -M 


(ni]+ni2)c^ ^ mi+ni2 


with 


fl-fl 

M = QtB + nil nij - r 

nil ~ ^2 


( 22 ) 


(23) 
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flYl “^7 flJl 

and Q = l + -• For the Chern-Simons Lagrangian given by eq.(3) 

ni] - m2 

above eq.(22) is the counterpart of eq.(46) in Ref.2 ;eq.(23) therefore defines the Runge- 
Lenz vector M of the Lagrangian (3) in the centre-of -momentum system. 

It is easy to see that the above derivation in (22) and (23) needs a second look for 
the equal mass mj = m 2 case .Indeed a similar discussion is also in order for the Runge- 
Lenz vector derived from eq.(45) in Dahl’s paper ^; note however that the remedy there 
is quite painless ,namely:begin with unequal masses mj ^ m 2 and subsequently derive the 
correct Lenz vector for the equal mass case, mj=m 2 =m, the result being a simple 
replacement of the factor p in eq.(46) there by m/2. But a corresponding effort here 
merits a separate discussion and is therefore relegated to the Appendix at the end of this 
paper. 

With Ml denoting the x-component of the 2-component vector M it is clear that 

d Ml . . 

^ ^ = QBj ;thus one expects the Poisson bracket [Mj ,H\ = -QBi , since M is a 


constant of motion,with H being given by eq. (9b). It is easy to check this using the 
following: 

{r,.fi] = ^^.{r,.f2] = ^^^.{r:.pj} = Sij.{f,.f2] = 0.{B.f,} = 0 = {B.f2] (24 

m; c tn 2 c 

On replacing r i by r 2 in the first pair of Poisson brackets in (24) one should also replace 
B 1 by B 2 . 

One can also use eq. (24) to evaluate the Poisson bracket of M with Lz the latter 
being given in eq. (9b). Using {b. ,L^} = -e^j Bj it is easy to verify that 
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( 25 ) 


with 812 = - 821 = 1 . We shall now evaluate the {Mj,M 2 } Poisson bracket ; by virtue 
of eq. (24) and the fact that [Bj, 82 ] = 0 it is easy to obtain the result 




with (nij - m 2 )cx = m^m 2 .While eq. (25) is expected on account of M being a vector, 


(26) does not match with that given by eq. (2) for the Runge-Lenz vector K and is 
therefore a novel feature of this paper. The fact that the calculation made here is exact 
and quite unlike the the approximate calculation to order 1/c of the Runge-Lenz vector by 
Dahl ^ reinforces our confidence in this assertion. 


The extra term eie2/7cc in eq. (26) is reminiscent of the observation by Witten and 
Ohve long ago that in supersymmetric theories with sohtons the usual supersymmetry 
algebra is modified to include topological quantum numbers as central charges. 
Interestingly ,the last factor on the right hand side of eq.(26) can be expressed as 

( I - 

^2 —1 —1 \^i I 

- -I- = id(rxl)^ where we define I = X]Y] -X 2 Y 2 with Yi = 1- — Y~ 



Quahtatively,eqs.(25) 


and 


(26) are important to this paper because 


we have derived them,besides eq.(24), as Poisson brackets here.This is not the case with 
the Poisson bracket relations given by eq.(4) in Ref.2; indeed Dahl provides only a 
quahtative understanding (see eqs.(56) - (58) in Ref.(2)) of the said Poisson bracket and 
is unable to derive eq.(4) in the 1/c order description in his paper.Needless to say so, 
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eq.(26) is the 2+1 dimensional counterpart of the 3 + 1 dimensional Poisson 


bracket given by eq.(4) in Ref.(2),namely 




2 

— H 


At k=l 


^ijk ^k 


(27) 


In conclusion ,two distinguishing features characterise the derivation of the Runge- 
Lenz vector (23) associated with the Chern-Simons Lagrangian (3) in this paper: a)The 
calculation here is exact and, (b) as an unexpected bonus the Poisson bracket in (26) 
contains the product eie 2 as a central charge. 
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When the masses m i and m 2 in the Lagrangian (3) are equal the Runge-Lenz vector M 
is worked out initially by assuming that nij m 2 ;taking the lim it of eq.( 22 ) when 

m 2 Mj then yields the desired form of M as explained below.Let us assume here that 

mj = m 2 +£ where 8 > O.For the second term in eq.(23) of the text one thus has 


lim-( fj - f2 ) = Urn < 

e^0£ e^O 


1 + 


B\ 


\l/2 ( 


c (m2 +£y 


1 + 




\l/2 


2 2 , 
C m2 J 


(Al) 


l^r 

^ jy2 

c m 2 


(A2) 


by a Maclaurin expansion of the first term to order 8 .Here 72 and yi (below) are given by 


eqs.(l 1) of the text. A similar effort on the first term in eq.(23) yields 


Urn 

e^0\ 


1- 


mifir2 -m2f2ri 


\B\ 

^72 2 2 

m2 C 


(A3) 


Thus when m 1 = m 2 = m ,we have instead of eq.(23) the Runge-Lenz vector 


M = 


B^y 


me 


2t 
— B-r 

V m 


A 


7(1-f )(mr-2tB) 


(A4) 


with f y =1 and 7 = 


1 + 


2 2 
m c 


With the Poisson brackets given by eqs.(24) in the 


text it is now easy to arrive at the relations [m i,L^] = -£ yMj and 

f 1 1 4 4 (^1^2 

[Mj,M2\ = -^rUf^-l\—+L, 


(A5) 


J 


Eq.(A5) is the equal mass counterpart of eq.(26) of the text for unequal masses.. 
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